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ABSTRACT 



Aims. Characterizing the gas and dust properties prior to and in the neighborhood of active intermediate- to high-mass star formation. 
Methods. Two Infrared Dark Clouds (IRDCs) - IRDC 19175-4 and IRDC 19175-5 - that are located in the vicinity of the luminous 
massive star-forming region IRAS 19175+1357, but that remain absorption features up to 70/jm wavelength, were observed with the 
Plateau de Bure Interferometer in the 3.23 mm dust continuum as well as the N 2 H + (l-0) and 13 CS(2-1) line emission. 
Results. While IRDC 19175-4 is clearly detected in the 3.23 mm continuum, the second source in the field, IRDC 19175-5, is only 
barely observable above the 3cr continuum detection threshold. However, the N2H + (l-0) observations reveal 17 separate sub-sources 
in the vicinity of the two IRDCs. Most of them exhibit low levels of turbulence (Av < 1 km s ), indicating that the fragmentation pro- 
cess in these cores may be dominated by the interplay of thermal pressure and gravity, but not so much by turbulence. Combining the 
small line widths with the non-detection up to 70yum and the absence of other signs of star formation activity, most of these 17 cores 
with masses between sub-solar to ~10M Q are likely still in a starless phase. The N9H + column density analysis indicates significant 
abundance variations between the cores. Furthermore, we find a large CS depletion factor of the order 100. Although the strongest 
line and continuum peak is close to virial equilibrium, its slightly broader line width compared to the other cores is consistent with 
it being in a contraction phase potentially at the verge of star formation. Based on the 3.23 mm upper limits, the other cores may be 
gravitationally stable or even transient structures. The relative peak velocities between neighboring cores are usually below 1 km s , 
and we do not identify streaming motions along the filamentary structures. Average densities are between 10 5 and 10 6 cirr 3 (one to 
two orders of magnitude larger than for example in the Pipe nebula) implying relatively small Jeans-lengths that are consistent with 
the observed core separations of the order 5000 AU. Environmental reasons potentially determining these values are discussed. 
Conclusions. These observations show that multiple low- to intermediate-mass low-turbulence starless cores can exist in the prox- 
imity of more turbulent active intermediate- to high-mass star-forming regions. While masses and levels of turbulence are consistent 
with low-mass starless core regions, other parameters like the densities or Jeans-lengths differ considerably. This may be due to 
environmental effects. The quest for high-mass starless cores prior to any star formation activity remains open. 
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1. Introduction 



Stars: individual: IRDC 19175-4, IRDC 19175-5 - Techniques: interferometric - Line: formation 



Although our understanding of massive star formation has 
improved tremendously over the last decade, most progress 
was made toward relatively evolved regions like High-Mass 
Protostellar Objects (HMPOs), Hot Molecular Cores (HMCs) 
or Ultracompact H n regions (UCHns), s e e for example severa l 
recent reviews like IBeuther et alJ (l2007l): iBonnell et al] (l2007h: 
ICesaroni et ail d2007l) : iHoare et all d2007l) : IZinnecker & Yorkd 
d2007l) . Toward even earlier evolutionary stages prior or at the 
onset of massive star formation, in recent years significant pro- 
gresses have been made in identifying large samples of Infrared 
Dark Clouds (IRDCs) via absorption studies i n the mid-infrared 
with the MSX ISO and Spitzer satellites ( e.g.,lEgan et all 19981 
iBacmann et ail l2000t ISimon et alJ 120061: iJackson et al] " l2008l) . 
Such IRDCs are expected to harbor genuine high-mass star- 
less cores (HMSCs) as well as high-mass cores with embedded 
low- t o intermediate-mass p r otostars in the early accretion phase 
(e.g., iForbrich et all l2004t iKlein et alJ 120051: ISridharan et ail 



20051: iRafhborne et al II2005L 120061: IBeuther & Steinac"kerll2007t 
Beut her & Sridharanl 120071) . While the samples contain more 
than 10000 IRDCs, it is interesting to note that from the pub- 
lished higher spatial resolution case studies (of the order 20) 
none is a real HMS C but all have ident i fied embedded accret- 
ing protostars (e.g.. IBeuther et al.l 12005]: IBeuther & Steinackerl 



2007; Birkmann et al. 2006; Ormel et al. 2005 
Rathborne et al.l 120051 l2008at IForbrich et al 



Pillai et al.l 2006: 
2009). Hence, a 



* The fits-files of the continuum and line images are available in 
electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http:// cdsweKu : strasbg.fr/cgi-bii Vq"cat?J/A+A/| 



genuine HMSC has not yet been unambiguously identified in 
the literature. 

To overcome this obvious missing link, here we report 
high-spatial-resolution observations and analysis of two excel- 
lent HMSC candidates in the dust continuum and the N2H + 
line emission to investigate their physical properties prior to 
any star formation activity. N2H + is known to be strong in 
low-mass starless core s (e.g., iBergin et"ai1 120021: iTafalla et"ail 
12001 IChen etaDl2007b as well as massive IRD Cs with already 
embedded, very young protostars (e.g., IBeuther et al.1 12005). 
Furthermore, it is a good chemical tracer of the cold early evo- 
lutionary phase because it does not deplete much onto grains 
at low temperatures but it starts getting destroyed when the core 
heats up releasing the CO from the dust grains (which was frozen 
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onto them during earlier cold phases) that subsequently react s 
with N2H + in the gas-phase chemistry (Ber gin & Langerlll997l) . 
While the dust continuum emission pinpoints the locations of the 
densest gas and dust components, the line emission will be used 
to determine the gas kinematics of the region. 

The two HMSC candidate sources are located in the 
vicinity of the High-Mass Protostellar Object (HMPO) 
IRAS 19175+1357. During our initial stud y of the dust c ontin- 
uum emission of a sample of HMPOs (Sridharan et al., 2002; 
Beut her et al.L 120021) . we serendipitously detected mm contin- 
uum peaks offset from the IRAS sources that turned out by com- 
parison with the MSX database to be associated with IRDCs 
dSridharan et al.Ll2.Q05l) . The average NH3 temperatures and line 
widths of these objects are significantly smaller than those from 
typical HMPOs and UCHns, hence putting them at an earlier 
evolutionary stage. However, follow-up SiO(2-l) observations 
revealed SiO line wing emission toward ~40% of the sample, 
indicating t hat star formation has starte d in at least that fraction 
of sources (Beu ther & SridharanTl2007h . 



grey: MSX at 8/im; contours: 1.2mm dust continuum 




R.A (J2000) 

Fig. 1. Large-scale emission in the environment of the sources 
IRDC 19175-4, IRDC 19175-5 and IRAS 19175+1357. The 
grey-scale shows the MSX 8 jum mid-infrared emission, 
and t he contours o utline the 1.2 mm continuum emission 
from IBeuther et alj d2002l) in steps of 20mJybeairT I (3<x 
rms~15 mJybeairT 1 ). The circle presents the PdBI primary 
beam size of ~ 54" encompassing the two target regions. A 
scale-bar and the IRAM 30 m beam size of the 1 .2 mm contin- 
uum observations are shown at the bottom left. 

To identify a candidate HMSC for this sample with the 
given datasets, we require the source to have (a) mm continuum 
emission, (b) no detectable SiO emission, (c) narrow NH? and 
H 13 CO + line width, (d) low NH3 temperatures and (e) no detec- 
tions in the Spitzer IRAC and MIPS 24 and 70 (im bands. Here, 
we select from the original sample two particularly promising 
HMSC candidates that fulfill these criteria and even fit within 
one PdBI 3 mm primary beam (Fig.Q]). 

The two sources IRDC 19175-4 and IRDC 19175-5 with 
a projected separation of ~22000AU are part of the 
larger-scale filamentary star-forming region associated with 
IRAS 19175+13 57 (Fig.Q) at an estimated distance of -1.1 kpc 
( Sridhara n et al.Ll2005l) . The luminosity of the main IRAS source 
in Fig. Q] is approximately 10 3 L G , lo wer than origina l ly esti - 
mated based on a wr ong distance in ISridharan et alj ([2002). 
Srid haran et al.l d2002l) found weak cm continuum emission in- 
dicative of ionized gas in the vicinity of the IRAS 19175+1357 
source. However, no maser emission or strong CO line wings 



were detected. The latter indicates either only weak or no out- 
flow activity anymore toward this main source, or the outflow 
could be oriented close to the plane of the sky. Since the vi sr 
of the main IRAS source and the two IRDCs are approximately 
the same (6.9 versus 7.7kms~ 1 ), we can safely assume that both 
regions are at about the same distance. Hence their projected 
separation of ~0.45 pc is a good approximation of their real sep- 
aration. Although we cannot proof that the IRAS source directly 
influences the two IRDCs, both regions should nevertheless stem 
from the same larger gas clump that originally fragmented into 
the already more evolved IRAS source and the younger IRDCs. 

Figure [2] prese nts an overlay of all IRAC and the M I PS 24 
and 70 //m bands dChurchwell et all 120091; ICarev et all l2009h 
with the dust continuum emission towards IRDC 19175-4 and 
IRDC 19175-5, and both regions remain dark even up to 70/rni. 
The 24/mi absorption feature toward IRDC 19175-5 is a bit 
more pronounced than toward IRDC 19175-4. We detected no 
SiO emission. Their NH3 single-dish line widths are ~1.5km/s, 
and their derived NH3 rotation temperatures are of the order 1 5 K 
dSridharan et all l2002t IBeuther & Sridharanl 120071) . All these 
characteristics make the two sources ideal HMSC candidates and 
hence well-suited targets to investigate the physical properties 
prior to the onset of star formation. 

Assuming optically thin emission from cold dust at 15 K and 
a dista nce of 1 . 1 kpc with a gas-to-dust mass ratio of 186 (fol- 
lowing D raine et"al]|2007l; lleii kins 2004), the integrated 1.2 mm 
emission fro m the two region IRD C 19175-4 and IRDC 19175-5 
of 620 mJy dBeuther et all 120021) corresponds to approximately 
87 M of gas mass, fairly equally divided between both sources. 
While these gas masses are not of the order 1OOOM like some 
better known massive star-forming regions, the gas reservoir is 
large enough to potentially form intermediate-mass B-type stars. 



2. Observations 

We observed IRDC 19175-4 and IRDC19175-5 with the Plateau 
de Bure Interferometer during two nights in August 2007 and 
April 2008 at 93 GHz in the C and D configurations covering 
projected baselines between approximately 13 and 180 m. The 
3 mm receivers were tuned to 92.835 GHz in the lower side- 
band covering the N2H + (l-0) and the 13 CS(2-1) lines as well as 
the 3.23 mm continuum emission. The phase noise was always 
lower than 30° and mostly lower than 20°. For continuum mea- 
surements we placed six 320 MHz correlator units in the band, 
the spectral lines were excluded in averaging the units to pro- 
duce the final continuum image. Temporal fluctuations of am- 
plitude and phase were calibrated with frequent observations of 
the quasars 1923-201 and 1827+062. The amplitude scale was 
derived from measurements of MWC349. We estimate the fi- 
nal flux accuracy to be correct to within ~ 15%. The phase 
reference center is R.A.[J2000] 19 A 19 m 50. 169 s and Dec.[J2000] 
14°01'13.75", and the velocity of rest vi sr is 7.7 km s . The syn- 
thesized beam of the continuum data with a robust weighting is 
3.5" x 2.8" (PA. 28° east of north) with a 3<x continuum rms 
of 0.16mJy beam -1 . The 3<x rms of the N2H + (l-0) data mea- 
sured from an emission-free channel with a spectral resolution of 
0.2 km s" 1 is 15 mJy beatrT 1 . The 3cr rms of the 13 CS(2-1) data 
measured from an emission-free channel with a spectral resolu- 
tion of 0.5 km s~' is lOmJy beatrr 1 . 

These mm observations are complemented with the mid- 
infrared SPITZER data from the GLIMPSE and MIPSGAL sur- 
veys of the Galactic plane using the IRAC camera centered at 
3.6, 4.5, 5.8 and 8.0/mi and the MIPS camera centered at 24 and 
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Fig. 2. Zoom into the central region of the two target sources. Like FigureQ] the contours show the 1 .2 mm dust continuum emission 
dBeuther et all 120021) . and the grey-scale presents the mid- to far-infrared emission as observed with the Spitzer Space Telescope. 
From top-left to bottom-right, the 6 panels show the Spitzer bands centered at 3.6, 4.5, 5.8, 8.0, 24 and 70/im, as marked in each 
panel. The dotted circles in each panel mark the only two Spitzer sources that were identified as classO/I candidates (see main text). 



70yum dWerner et ail l2004t iFazio et aTl l2004t iBeniamin et all density sensitivities of 0.9 M Q and 1.7xl0 23 cm 2 . Comparing 
120031: iRieke et all 12004 ICarev et all 120051 



3. Results 

3. 1. Millimeter continuum emission and class 0/1 sources 
identified via the Spitzer data 

Figure |3]presents an overlay of th e 1 .2 mm dust contin uum data 
from the IRAM 30 m telescope dBeuther et all |2002|) with the 
new high-spatial-resolution 3.23 mm data obtained now with the 
PdBI. While the north-eastern peak (IRDC 19175-4) is clearly 
detected also in the new 3.23 mm data, the south-western peak 
(IRDC 19175-5) is a marginal detection at less than 5<x with 
0.26mJybeam _1 . The peak and integrated 3.23 mm fluxes of 
IRDC 191 75-4 are l.OgmJy bea m" 1 and 1.8 mJy, resp ectively. 
Following Iffildebrandl (1 19831) and lBeuther et all (120021) . assum- 
ing optically thin d ust continuum emission with a g as-to-dust 
mass ratio of 186 ( Dra ine et all 120071: iJenkinsl |2004|) . a dust 
opacity index ( 3 = 2 (corresponding to a typical dust distribution 
as outlined in iMathis et al.lll977t lOssenkopf & Henn ing 1994) 
and a temperature of 15 K dSridharan et all | 20051) . we derive a 
gas mass and column density of ~1OM and ~1.8xl0 24 cm~ 2 , 
respectively . The latter correspon ds to a visual extinction A v of 
-1900 mag dFrerking et all [T982). Our 3cr continuum sensitivity 
corresponds with the given parameters to 3<x mass and column 



the mass derived from the PdBI data of IRDC 19175-4 of 
~1OM with the single-dish estimated mass of 87 M , approxi- 
mately 90% of the flux has been filtered out by the interferometer 
observations. 



Although the mm continuum peaks are infrared dark up 
to 70jum emission, Figure [2] shows several Spitzer sources 
within the field of view. Empl oying the available IRA C 
GLIMPS E point source c atalog dChurchwell et alj 120091 & 
http://irsa.ipac.caltech.edu/) and adopting the Spitzer color cri- 
teria by lAllen et aTl d2004h . we tried to identify possible class 
0/1 and class II protostellar candidates around the discussed 
IRDCs (see also Fang et al., subm.). Interestingly, we do not 
find any class II source and only 2 class 0/1 candidates close to 
IRDC 19175-4 and IRDC 19175-5 (see dotted circles in Fig.©. 
While these two class 0/1 candidates are clearly within the vicin- 
ity of the IRDCs, almost all of the remaining sources in the field 
of view are either more evolved and hence do not fit the Spitzer 
color criteria, or they are entirely unrelated fore- and/or back- 
ground sources. Identifying only two young protostellar candi- 
dates with locations not toward the mm continuum peaks further 
stresses the youth of the target regions. 
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Fig. 5 . Channel map of the isolated N 2 H + (l-0) (Fl, F=0,l-»1,2) hyperfine emission component at 93176.2650MHz (ICaselli et all 
1995) which is blue-shifted by -8.0kms _1 with respect to the central component (see Fig. [7J»- The velocities as marked in the top- 
right corners of each panel are shifted to the velocity of rest. The grey-scale and contours are done in 3<x steps of 15 mJy beam -1 
(full and dashed contours positive and negative, respectively). The two stars mark the main two mm continuum sources (Figs. [TJ 
& [3), the numbers in the top-right corners show the central velocities of each channel, and the synthesized beam is plotted at the 
bottom-left of each panel. The offsets are given with respect to the phase center. 



grey: 1.2irim continuum (MAMBO) 
contours: 3.23mm continuum (PdBI) 
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Fig. 3. Millimeter continuum emission toward IRDC 19175- 
4 and IRDC19175-5. The grey-scale with dotted contours 
presents the single-dish 1 .2 mm continuum observations with the 
MAMBO array mounted on the IRAM 30 m telesc ope contoured 
at the 3cr levels of 15mJybeam _1 (FigJTJ and iBeuther et all 
|2002|) . and the thick contours show the new 3.23 mm data 
from the PdBI observed for this project with 3cr contours of 
0.16 mJy beam - 1 . Full and dashed contours present emission and 
negative features due to the insufficient uv-coverage, respec- 
tively. In the bottom-left corner, the beam of the 1 .2 mm single- 
dish observations is shown as the large circle, and the synthe- 
sized PdBI beam as smaller ellipse. 



3.2. N 2 H + (1-0) emission 

Figures [4] [5] and [6] present the N2H + (l-0) data as integrated 
intensity maps, a channel map and the deduced 1st and 2nd 



moment maps (intensity-weighted peak velocities and line 
widths). For the channel and moment maps, we use the isolated 
N2H + (l-0) (Fl, F =0, 1 — » 1 ,2) hyperfin e emission component at 
93176.2650MHz dCaselli et all 1 19951) which is blue-shifted by 
-8.0km s with respect to the central component because it is a 
spectrally isolated line without further line-blending and mostly 
optically thin. In strong contrast to the relatively simple dust con- 
tinuum structure, the N2H + data now reveal a large amount of 
separate emission peaks, 17 in number. Although the general 
structure of the two integrated maps of the main and satellite 
hyperfine structure lines in Figure [4] is very similar, small dis- 
crepancies arise because of different sensitivities and opacities. 
While N2H + peaks 10, 11, 14, 15 and 17 are mainly detected in 
the main hyperfine structure line most likely because it can trace 
lower column density structures, other peaks are offset between 
the two integrated images. For example the N2H + peaks 7 and 8 
are clearly separated in the satellite hyperfine structure line map, 
whereas it appears like a single peak in the main hyperfine line. 
This difference likely arises from the higher optical depth of the 
main line. For our source-labeling we refer to the positions from 
the isolated satellite hyperfine structure line where they were 
clearly detected. For the remaining peaks, the positions are ex- 
tracted from the main hyperfine line. While the strongest N2H + 
peak number 1 is associated with the strongest mm continuum 
peak, the second mm peak has no directly associated N2H + fea- 
ture, rather the N2H + peaks 4 and 5 are lying at the edge of the 
mm continuum emission. FigureQpresents a few example N2H + 
spectra, and Table QJ shows the corresponding spectral parame- 
ters (peak line temperature r pea k, peak velocity v, line width Av 
and total optical depth t of all hyperfine components) we de- 
rived from fitting the complete N2H + (l-0) hyperfine structure 
within CLASS (part of the GILDAS software packageQ). While 
the N2H + peak 16 is spatially associated with the western Spitzer 
class 0/1 source (Figs.|2]&|6]l, no other N2H + peak has a Spitzer- 
identifiable young protostellar counterpart. Except of peak 16, 
all N2H + peaks are good candidates for being genuine starless 
cores. 



http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 6. Moment maps of the N2H + (l-0) emission. The grey-scale in the left and right panels present the 1st and 2nd moments, 
respectively (intensity-weighted peak velocities and line widths). The contours in both panels show the Oth moment or integrated 
intensity maps which were produced from the channel map in Fig.[5]via clipping all data below 5cr corresponding to 25 mJy beam" 1 
for a channel width of 0.2km s -1 . The contour levels are done in steps of lOmJybeam" 1 kms" 1 . The maps were produced from 
the isolated N2H + (l-0) (Fl, F=0,1— >1,2) hyperfine emission component (-8.0kms _1 offset from the central component) with an 
integration regime corresponding to 6.3 to 8.9 km s" 1 (vi sr = 7.7 km s ). The stars mark the positions of the main mm continuum 
peaks (Figs. [T] & O, and the synthesized beam and scale-bar are shown at the bottom-left of each panel. The numbers in the left 
panel mark the 17 N2H + (l-0) emission peaks discussed in the main text and Table Q] The dotted circle associated with the N2FF 
peak 16 marks the position of the western Spitzer class 0/1 source. The offsets are given with respect to the phase center. 



The sub-structure of the cores exhibits an elongated and fila- 
mentary structure. The filament associated with the N2FF peaks 
6 to 10 has a linear extend of ~25000 AU with a length-to-width 
ratio exceeding 5. The a verage separation betw een the five cores 
is -5000 AU. Following lStahler & Pallal (120051) we can calculate 
the Jeans-length 



^ = (ff!)" 2 = „. 19pc (_L)" 2 (^_)-" J 

\Gp) F \10K/ \10 4 cm" 3 / 

with the sound speed a, the gravitational constant G and the 
density p. What density do we expect for the cores? The critical 
density of the N2H + (l-0) transition at the given low temperature 
is ~ 1.5 x 10 5 cm" 3 . In addition to this, we can calculate aver- 
age densities for each core converting the N2H + column densi- 
ties (see section [3~3l to H2 column densitie s assuming a typica l 
N 2 H + to H 2 abundance ratio of 3 x 10" 10 dCaselli et al.ll2002al) . 
and furthermore assuming in spherical geometry that the 3rd di- 
mension along the line of sight corresponds to the synthesized 
beam for the compact cores. Table Q] lists the derived densities 
for all cores estimated this way. The range of average densities 
varies by about one order of magnitude between 2 x 10 5 and 
1 .9 x 10 6 cm" 3 . For the cores 6 to 10 along the filament, the aver- 
age density is about ~ 6 x 10 5 cm" 3 , implying that the peak den- 
sities should be even higher. Hence, we can estimate the Jeans- 
length at densities of 5 x 10 5 and 10 6 cm" 3 with a typical temper- 
ature of 15 K. This brackets the corresponding Jeans-lengths ap- 
proximately between ~4800 and ~6800 AU above which the fil- 
ament should become unstable. Interestingly, the measured core 
separation corresponds well to this Jeans-length indicating that 
the cores could be stable against further fragmentation. 



While the broadest line (still only 1.8kms"' full width half 
maximum) is associated with the main mm continuum peak, 
all other peak positions exhibit line widths close to or be- 
low lkms" 1 . The N 2 H + (l-0) thermal line width at 15K is 
-O.lSkms" 1 . Since the spectral resolution of our final N2H + (1- 
0) data cube (0.2 kms" 1 ) is insufficient to actually resolve the 
thermal line width, it is difficult to accurately constrain whether 
the narrowest measured line widths (Table [TJ close to our spec- 
tral resolution (less than 3 spectral resolution elements across 
the line width Av) are not rather upper limits to real line widths 
closer to thermal. 

The relative peak velocities between close-by neighboring 
sub-cores is usually below 1 kms" 1 , and we do not identify any 
clear velocity gradient or streaming motions along the filamen- 
tary structures associated with the N2I-F peak positions (see 
Figure [6] and Tabl^B- 

3.3. N 2 H + column densities and abundances 

Using the line intensities and optical depths from our N2H + (1- 
0) fitting, we can estimate the N2H + column densities for all 
cores (e.g. ICaselli et al.ll2002al) . The resulting values between 
0.3 x 10 13 and 6.6 x 10 13 cm" 2 are also listed in Tabled] With 
a typic al N 2 H + /H 2 abundance ratio of 3 x 10" 10 (ICaselli et all 
l2002al) . this corresponds to an H2 column density regime be- 
tween roughly 2 x 10 22 and 2 x lO^cm" 1 (or v isual extinc- 
tions A v between 20 and 200. iFrerking et al.|[l982b . largely be- 
low our H2 3<x column density limit of 1.7 x 10 23 cm" 2 derived 
from the 3.23 mm dust continuum emission (see section [3~TT >. 
However, taking our derived H2 and N2H + column densities at 
face value, we derive an order of magnitude lower abundance 
of ~ 3.7 x 10" 11 for the main 3.23 mm continuum and N2H + 
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Table 1. Parameters of N2H + (l-0) emission peaks 



# 


Offset" 


^peak 


V 


Av 


T 


AW 


M*. 

vir 


av. p 






K 


km/s 


km/s 




xl0 13 cnr 2 


M G 


10 5 cirr 3 


1 


6.2/8.6 


2.1 


7.7 


1.8 


11.9 


6.6 ±0.3 


7.5-11.4 


19 


2 


12.1/17.8 


0.7 


6.9 


0.8 


12.7 


1.9 ±0.3 


0.6-0.9 


12 


3, 


0.7/4.3 


1.5 


7.3 


0.5 


7.0 


0.9 ±0.1 


0.3-0.4 


5.9 


3 2 


0.7/4.3 


0.4 


8.3 


0.8 


16.3 


2.1 ±0.5 


0.6-0.9 


14 


4 


-10.1/-9.0 


2.3 


8.1 


1.1 


6.5 


2.3 ± 0.2 


1.3-1.9 


15 


4a 


-9.9/-5.8 


2.3 


7.7 


0.7 


7.8 


1.7 ±0.1 


0.5-0.7 


11 


5 


-8.6/-4.5 


2.2 


7.7 


0.8 


7.1 


1.9 ±0.1 


0.7-1.1 


12 


6 


-5.7/0.7 


1.3 


8.2 


0.7 


12.6 


2.0 ±0.1 


0.5-0.7 


13 


7 


-3.6/5.3 


1.3 


8.1 


0.6 


5.1 


0.7 ±0.1 


0.3-0.5 


4.6 


8 


-2.9/8.2 


4.3 


8.0 


0.7 


0.1 


0.8 ±0.1 


0.5-0.7 


5.2 


9 


-1.0/13.0 


0.9 


8.2 


0.5 


8.3 


0.9 ±0.2 


0.3-0.4 


5.9 


10 


-0.1/17.8 


2.2 


8.1 


0.5 


2.0 


0.4 ±0.1 


0.3-0.5 


2.6 


Hi 


4.6/23.7 


1.4 


7.6 


0.5 


1.1 


0.3 ±0.1 


0.3-0.5 


2.0 


H 2 


4.6/23.7 


0.6 


8.6 


0.8 


4.9 


0.7 ±0.1 


0.6-0.9 


4.6 


12 


-14.2/- 11.1 


2.0 


7.7 


1.2 


5.7 


2.0 ±0.2 


1.5-2.2 


13 


13 


-S.4/-7.4 


1.5 


7.8 


0.7 


10.4 


1.8 ±0.1 


0.5-0.7 


12 


14, 


-0.6/-5.2 


0.5 


7.8 


1.0 


10.7 


1.9 ±0.2 


1.0-1.5 


12 


14, 


-0.6/-5.2 


2.6 


8.8 


0.7 


0.1 


0.3 ±0.1 


0.5-0.7 


2.0 


15 


-17.6/-4.9 


2.6 


8.9 


0.5 


0.9 


0.3 ± 0.3 


0.3-0.4 


2.0 


16 


-24.5/3.3 


2.8 


6.8 


1.0 


3.7 


1.4 ±0.1 


0.9-1.4 


9.2 


17 


-20.7/5.7 


3.1 


6.8 


1.0 


0.1 


0.6 ± 0.2 


1.0-1.6 


3.9 



The subscripts associated with sources 3, 1 1 and 14 correspond to two different velocity components required for the spectral fits. 
" The offsets are given with res pect to the phase center R.A.[J20001 19*19"*50.169 s and Dec.[J2000] 14°01' 13.75". 
b The virial masses were calculated following MacL aren et alj ( [19881) for 1/r and 1/r 2 density distributions, resulting in the ranges of values. All 
source radii were approximated with half of the average NiH + synthesized beam of ~ 3.1". Only for peak one we used a source radius of 3.5". 



peak 1 associated with IRDC 19175-4. Assuming that this abun- 
dance could hold for all N2H + (l-0) cores, the range of cor- 
responding H2 column densities would be an order of mag- 
nitude higher, roughly between 2 x 10 23 and 2 x 10 24 cm~ 2 . 
These values nearly all were above our H2 3<x column density 
limit of 1.7 x 10 23 crrT 2 derived from the 3.23 mm continuum 
emission. Although all these estimates are prone to many un- 
certainties, they indicate that even during such early evolution- 
ary stages the gas phase abundances are not uniform through- 
out the core. The source with the broadest line width, and hence 
likely the most evolved of the cores (peak 1 in IRDC 19175-4), 
shows the lowest a bundance compared to t he other less turbu- 
lent cores. While Tafal la & Santiago! ( 120041) argued in the case 
of the low-mass starless core LI 57 IE that particular low N2H + 
abundanc es could also be due to very young evolutionary stages 
(see also [Aikawa et al. 2003), the increased line width toward 
our main core makes this an unlikely scenario. Rather in con- 
trast, increased densities in more evolved stages of core collapse 
(like those inferred here, see TableQ]) increase t he recombination 
rate a nd hence lower the N2H + abundance (e.g, Bergin & Langer 
Il997[) . This appears the more likely scenario for the cores we ob- 
serve here. 

3.4. n CS(2-1) emission 

In addition to the N2H + (l-0) emission, our spectral setup also 
covered the 13 CS(2-1) line. However, although detected, the 
13 CS(2-1) emission is significantly weaker than N2H + (l-0). 
Figure [8] presents an overlay of the integrated 13 CS(2-1) with 
the N2H + (l-0) emission. It is interesting to note that while the 
main mm and N2H + peak is also detected in the 13 CS data, of the 
remaining N2H + cores, only peak 6 is associated with detectable 
13 CS emission above the 3cr level. This already indicates pe- 
culiar low abundances compared to more evolved regions, but 
likely also varying abundances between the different cores. 



Toward the main continuum and line peak, the integrated 
13 CS(2-1) emission of 0.042 Jy km s _1 corresponds in the 
Rayleigh-Jeans limit to 0.52 K km s . How does that compare 
to total single-dish fluxes? iBeuther et alj (2002) observed the 
CS(2-1) line toward that region with the IRAM30m telescope, 
and the measured integrated intensity of 3.3Kkms _1 corre- 
sponds to a flux of 17.6Jykms _1 . Assuming a typical line in- 
tensity ratio between 12 CS and 13 CS of ~7 we get an integrated 
13 CS intensity of 2.45 Jy km s _1 (toward the hot core NGC6334I 
and the much younger sibling NGC6334I(N), the measured ra- 
tios are 8 and 6.5, Walsh et al. in prep.). Therefore, the missing 
flux we loose by the interferometer observations amounts to ap- 
proximately 80%, similar to the ~9Q% for the mm continuum 
data (see section [XT). 



Assuming optically thin emission, we can use the CS(2- 
1) data to estimate 13 CS column densi t ies of the compact gas 
component following R ohlfs & Wilson! d2006l) at the given low 
temperatures of ~15 K. Employing further a 12 CS/ 13 CS ratio of 
50, we can estimate the main isotopologue's 12 CS column den- 
sities and derive the CS abundance by comparison with the H2 
column densities calculated in section 13.11 This results in a CS 
column density of ~ 1 x 10 14 ctrT 2 and a CS abundance with 
respect to H2 of ~ 5 x 10 _1 1 . Such very low values are consistent 
with chemical models that predict CS colum n densities of that 
order at the beginning of core evolution (e.g.. lNomura & M illar 
|2004|) . The estimated abundance is about 2 orders of magni- 
tude lower than observed fo r example toward Orion-KL (e.g., 
Ivan Dishoeck & Bla ke 1998), indicating CS depleting factors of 
~ 100 in very young cores. For comparison, IZhangetalJ ([2009) 
recently found similarly high CO depletion factors (between 100 
and 1000) toward the IRDC G28.32. 
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Fig. 4. The thick contours present the integrated emission of the 
isolated N2H + (l-0) (Fl, F= 0, 1 — » 1 ,2) hyperfine emission com- 
ponent at 93176.2650MHz dCaseUi et all 1 1995b which is blue- 
shifted by -8.0 km s _1 with respect to the central component (see 
Fig. |7]i. The contour levels are in steps of lOmJybeanT 1 kms -1 . 
The grey-scale with dotted contours shows the integrated emis- 
sion of the central few hyperfine components around the 
vi sr with contour levels of 20, 40 and then continuing in 
40 mJy beanT 1 kms -1 steps. To produce the integrated intensity 
(or 0th moment) maps, data were clipped below 5<x correspond- 
ing to 25 mJy beatrT 1 for a channel width of 0.2kms _I . The 
corresponding integration regimes are 6.3 to 8.9kms _I for the 
isolated hyperfine structure line and 5.5 to 10.7 km s 1 for the 
blended central line, respectively. The stars mark the positions 
of the main mm continuum peaks (Figs. Q] & |3), and the syn- 
thesized beam and scale-bar are shown at the bottom-left. The 
numbers mark the 17 N2H + (l-0) emission peaks discussed in 
the main text and Table Q] The dotted circle associated with the 
N2H + peak 16 marks the position of the western Spitzer class 0/1 
source. The offsets are given with respect to the phase center. 



3.5. Virial masses 

Following Mac Laren et ail dl988l) . we calculated the virial 
masses M v ; r = £2 x R x Av 2 with the constant £2 between 190 
and 126 for 1/r or l/r 2 density distributions, the core radius R 
(in units of parsecs) and the N2H + (l-0) line width Av. For the 
radius we assumed almost all sources to be unresolved with an 
average N2H + (l-0) synthesized beam of ~ 3.1", corresponding 
to the core diameter. Only peak 1 is resolved and we used a ra- 
dius of 3.5". TableQ]lists the range of virial masses for the given 
line widths and assumed density distributions. 

The virial mass derived for the main core (between 7.5 and 
11.4M Q ) corresponds well to the gas mass derived above from 
the mm continuum emission (see section [3~TT i. This indicates that 
the core may be in virial equilibrium. Nevertheless, since this 
core also exhibits the broadest line width, it is possible that the 
core may be in a contract ing phase prior to active star formation 
(e.g. JCaselli et al.ll2002bl) . We doubt that the line broadening is 
due to an early outflow si nce we did not detect any SiO emission 
in previous observations dBeuther & Sridharanll2007l) . For most 
other N2H + cores, we cannot draw conclusions regarding their 




5 , )Q 15 

velocity (km/sj 

Fig. 7. Example N2H + (l-0) spectra extracted toward the posi- 
tions labeled in each panel (see Figure [6] and Table [T). The full 
lines show the fits with the parameters listed in Table Q] The y- 
axis shows the brightness temperature. 



virial status since most of the virial masses are at or below our 
3cr mass sensitivity of ~0.9 M derived from the continuum data. 
Only the virial masses of the N2H + cores 4, 10 and 12 are above 
this 3cr sensitivity limit, indicating that at least these cores are 
stable against gravitational collapse or may even be transient, 
again expanding structures. 
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Fig. 8. The grey-scale shows the 0th moment or integrated inten- 
sity N 2 H + (l-0) map of the isolated N 2 H + (l-0) (Fl, F=0, 1->1,2) 
hyperfine emission component as in Figure [6] The thick con- 
tours present the 13 CS(2-1) emission integrated from 6.9 to 
9.5 km s _1 (0th moment). The contouring is done in 3<x steps of 
5 mJybeam" 1 kms" 1 . The stars mark the positions of the main 
mm continuum peaks (Figs. [T]&[3]), and the synthesized beam 
is shown at the bottom-left. The offsets are given with respect to 
the phase center. 



4. Comparison with other star-forming regions 

How do the results of these two IRDCs compare with more typ- 
ical regio ns of low-mass s tarless cores, for e xample the Pipe 
nebula dAlves et al.L 120071: lLada et al.L l2008t IRathborne etafl 
2008b)? While some parameters like the masses and the ma- 
jority of cores being starless are similar between the Pipe and 
the IRDCs discussed here, other parameters are not. For exam- 
ple, the mean density of the Pipe cores is 7 x 10 3 cm 4 , approx- 
imately 2 orders of magnitude lower than what we estimate for 
the IRDC 19175 cores. This density difference results also in dif- 
ferent Jeans-length of the order 0.2 pc for the Pipe and ~6000 AU 
for the cores discussed here. Interestingly, like in our case where 
the observed core separation corresponds well to the Jeans- 
length, the same is the case for the Pipe nebula with a much 
wider average core separation of ~0.26pc (Joao Alves, private 
communication). This supports the long-standing argument that 
the relatively simple Jeans-analysis can be considered a good ap- 
proximation for the expected fragmentation scales. But why are 
the average densities so different? For a more conclusive answer 
to this question, a larger statistical sample would be required, 
however, the environments of the two regions already indicate a 
possible reason. While the Pipe nebula is relatively isolated, the 
IRDCs discussed in this paper are close to the intermediate- to 
high-mass star-forming region IRAS 19175+1357. Therefore, it 
is likely that the environment - either direct influence from the 
more evolved IRAS source and/or the density and column den- 
sity structure of the original gas clump - influences the density 
structure of the nearby cores and this way also determines the 
Jeans- and fragmentation length of the neighboring next genera- 
tion of star formation. 



The observed line width in typical low-mass starless cores 
is usually found to be close to thermal where the non-thermal 
line width contribution rarely exceeds the thermal part (e.g . , 
lAndre et all 120071: IRathborne etafl l2008bt iFoster et alj [2009). 
The situation is slightly different in starless cores associated with 
high-mass star-forming regio ns where on average broader line 
widths are found. For example. lFontani e t al. (2008) found N 2 H + 
line width around 0.7kms _1 in the intermediate- t o high-mass 
star-fo rmin g region IRAS 053 45+3157, whereas I Wang et aU 
(2008D and IZhang etaT] d2009t) found NH 3 line widths around 
1.7kms _1 towards the youngest massive core associated with 
G28.34+0.06. 

The average line width found in the 17 cores of IRDC 19175- 
4 and IRDC-19175-5 is ~0.8kms _1 with a spread between 0.5 
and 1.8 km s -1 . As discussed in section [3~2l the lower end of 
that regime is less accurately determined since our velocity res- 
olution is only 0.2 km s" 1 , bar ely resolving the narrowest lines. 
Similar to Font ani et alj d2008l) . our measured line widths exceed 
those of the low-mass cores, however, they are still narrower 
than observed toward more evolved massive star-forming re- 
gions (e.g., high-mass pro toste liar objects, HMPOs, orultracom- 
pact Hn regions, UCHns, iBeuther et alj|2002t [Churchwell et all 



1990; Cesaron i et al.ll 1991 ). This indicates a lower level of tur- 
bulent motions at the onset of intermediate-mass star formation 
compared to its successive evolution. Furthermore, the many 
N 2 H sub-sources show that such intermediate-mass cores can 
fragment significantly and potentially form small groups or clus- 
ters. 

Maybe a bit surprising, the line widths found towa rd the very 
young clump in G28. 34+0.06 by IWang et all d2008l) are inter- 
mediate between the narrower lines we and[Fontani et al. (2008) 
observed, whereas they are still less broad than than those of 
typical HMPOs or UCHn regions. It is interesting to compare a 
few o ther quantities of this particular region PI in G28.34+0.06 
(see IZhanget al.1 12009): The average separation between their 
5 sub-cores is 0.19pc, similar to the Pipe and different to the 
case discussed here. In contrast to that, the average densities in 
G28.34+0.06 range between 10 6 and 10 7 cm~ 3 , comparable to 
the values found in IRDC 19175-4 and IRDC-19175-5. As out- 
lined above, for the Pipe nebula, densities, Jeans-length and ob- 
served core separation match rather well. This is not the case 
for the PI region in G28.34+0.06. Their densities imply smaller 
Jeans len gth and lower Jean s-masses than observed. As dis- 
cussed by IZhang etaT] d2009), in the case of G28. 34+0.06, not 
only thermal pressure and gravity seem to control the fragmen- 
tation of the initial gas clump, but other processes are required. 
The broader observed line width there indicates a significant tur- 
bulent contribution which may help to stabilize the cloud cores 
against further fragmentation at that evolutionary stage. 

The situation for IRDC 19175-4 and IRDC-19175-5 is dif- 
ferent to G28.34+0.06. Since here densities, Jeans-lengths and 
observed core separation match well, turbulence seems to play 
a less important role for the fragmentation processes. The cloud 
sub-structure can well be dominated by the interplay of mainly 
thermal pressure and gravity. 

5. Discussion and Conclusions 

Combining Spitzer IRAC and MIPS data from 3.6 to 70jum 
with mm continuum and N 2 H + observations from the Plateau 
de Bure Interferometer and the IRAM 30 m telescope, we 
identify two very young gas and dust cores of intermediate 
mass (~87M G together for IRDC 19175-4 and IRDC 19175- 
5) in the vicinity of the intermediate- to high-mass star- 
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forming region IRAS 19175+1357. The high-spatial-resolution 
PdBI 3.23 mm continuum data clearly resolve IRDC 19175-4. 
However, IRDC 19175-5 is only barely detected above the 3<x 
level. More importantly, the N2H + (l-0) line reveals a com- 
pletely different structure with 17 separate emission sources 
within the vicinity of IRDC 19175-4 and IRDC 19175-5. Fitting 
the full N2H + (l-0) hyperfine structure, we are able to derive the 
line parameters (intensity, peak velocity, line width and optical 
depth) for all cores, and from that the N2H + column densities 
and virial masses. 

Of particular interest are the measured line width of the 
N2H + cores. Only the most massive core IRDC 19175-4 with a 
mass of ~1OM exhibits a line width of -l.Skms" 1 whereas 
all other N2H + cores have significantly narrower line widths 
below or around lkms -1 . While the measured line widths 
between 0.5 and 1.8kms _1 are still larger than the thermal 
value of ~0.15kms _1 (at 15 K), the line widths are narrow 
comp ared to more evolved high-mass star-forming regions 
(e ^lChurchwell et alil990t lHatchell et al]fl998l: iBeuther et alJ 
2002) and hence indicate a low level of internal turbulence. 
Matching densities, Jeans-lengths and observed core separation 
indicate that the fragmentation process may be dominated by the 
interplay of mainly thermal pressure and gravity, but not so much 
by turbulence. These low-turbulence cores are in the vicinity of 
a high-mass star-forming region with approximately 9 x 10 4 AU 
separation from IRAS 19175+1357. As outlined in section [1] 
although we cannot proof direct interaction between the more 
evolved IRAS source and the younger IRDCs, both should stem 
from the same gas clump that originally fragmented and pro- 
duced (massive) star-forming regions of different evolutionary 
stages. 

A comparison of the N2H + column densities with the cor- 
responding H.2 values (and upper limits) derived from the dust 
continuum emission, indicates significant abundance differences 
between the N2KF cores. While most of the lower-mass cores 
are consistent with typical N2KF abundances of the order 10~ 10 
with respect to H2, the main core IRDC 19175-4 has lower abun- 
dances of the order 10 . Since this source also shows the 
broadest line width, we interprete this lowered abundances as 
likely due to higher densities from the contracting core which in- 
creases the recombination rate and hence lowers the abundances 
of the ions. Furthermore, the weak 13 CS(2-1) detection allows 
us to estimate CS column densities and abundances, and we find 
CS depletion factors of the order 100 compared to more evolved 
regions like Orion-KL. 

The virial analysis and line width measurements indicate that 
the main core IRDC 19175-4, although being close to virial equi- 
librium, may already be in a phase of core contraction and hence 
at the verge of star formation. For the remaining cores, the sta- 
tus can be different. Since we have only upper limits for the gas 
masses from the dust continuum emission, a real virial analysis 
is difficult, however, these upper limits are consistent with virial 
stable cores. Similarly, the separation between the different sub- 
cores corresponds approximately to the Jeans-length at 15 K and 
a density of 10 6 ctrT 3 . Hence, the cores can be stable against 
further fragmentation. All this indicates that many of them may 
never enter star formation activity but could rather be transient or 
stable cores in the vicinity of other massive star-forming regions. 

In summary, here we report the characterization of 17 low- 
to intermediate mass N2H + cores with very low levels of tur- 
bulence in the vicinity of a high-mass star-forming region. The 
non-detection up to 70 //m wavelength, the low levels of turbu- 
lence as well as no other signs of star formation in previous 
studies indicates that we are dealing with starless cores. While 



some of them maybe contracting and may start active star forma- 
tion in the future, many of them may be stable against collapse 
and hence potentially even only transient structures. The average 
densities are approximately 2 orders of magnitude higher than 
in the Pipe nebula which is also manifested in smaller core sep- 
arations. We note that these rather quiescent cores are spatially 
associated with more turbulent intermediate- to high-mass star- 
forming cores, implying that low and high turbulence regions 
can co-exist in relatively close proximity. 
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